EVIDENCE FOR NEARBY SUPERNOVA EXPLOSIONS 
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Supernova explosions are one of the most energetic — and potentially lethal — phenomena in the 
Universe. Scientists have speculated for decades about the possible consequences for life on Earth 
of a nearby supernova, but plausible candidates for such an event were lacking. Here we show that 
the Scorpius-Centaurus OB association, a group of young stars currently located at ~ 130 parsecs 
from the Sun, has generated 20 SN explosions during the last 11 Myr, some of them probably as 
close as 40 pc to our planet. We find that the deposition on Earth of ^''Fe atoms produced by 
these explosions can explain the recent measurements of an excess of this isotope in deep ocean 
crust samples. We propose that ~ 2 Myr ago, one of the SNe exploded close enough to Earth to 
seriously damage the ozone layer, provoking or contributing to the Pliocene-Pleistocene boundary 
marine extinction. 

PACS numbers: 97.60.Bw, 26.30.+k, 91.50.-r, 98. 38. Am, 87.50Gi, 87.23Kg 



It has been proposed that the Local Bubble, a 150 pc 
hot (T « lO^K), low-density gas cavity which surrounds 
the solar system, was formed by several SN explosions 
during the last « 10 Myr [Q. The paucity of SNe in 
the Galaxy makes very unlikely that several isolated SN 
explosions would happen in short succession within such 
a small region, but about 20% of all SNe originate in OB 
star associations, and are therefore strongly clustered in 
time and space. By tracing back in time the positions 
of all nearby OB associations, it is possible to show that 
the Sco-Cen association was the only one able to produce 
SNe in the right numbers and places to generate the Local 
Bubble @. 

This association can be divided into three subgroups 
1^: Lower Centaurus Crux (LCC), Upper Centaurus Lu- 
pus (UCL) and Upper Scorpius (US). Using detailed age 
Q and membership information [|| , it is possible to com- 
pare the current numbers of early/late OB stars in each 
subgroup with the predictions of synthetic star formation 
models. The agreement with these models is good and 
based on them it can be estimated that each subgroup 
started producing SNe at 3-5 Myr after their formation, 
i.e. 7 — 8, 10—11 and 2 — 3 Myr ago respectively, with 
an expected constant rate of ^ 1 SN Myr~^ (see Ref. ||] 
for details). 

The current distances to the Sco-Cen subgroups can 
be reliably calculated using trigonometric parallaxes 
and then traced back in time taking into account the 
motion of both the Sco-Cen association and the Sun with 
respect to the local standard of rest which rotates with 
the Galaxy, as described in Ref. |^ . Fig 1 . shows how the 



distance from Earth of each subgroup has evolved during 
the epoch in which they were actively producing SNe. At 
its closest, about 2 — 3 Myr ago, the center of LCC was 
at ^ 100 pc from the Solar System; the spatial extent 
of these groups can be approximated by a Gaussian with 
cr = 25 — 30 pc, what means that SNe from LCC could 
have exploded as close as ~ 40 pc from Earth {2a lower 
limit). 

The explosion of a nearby SN can be detected by iso- 
tope anomalies in the geological record caused by the 
deposition of SN debris Recently, Knie et al. |^ 
measured a significant excess of ^°Fe atoms in two lay- 
ers of deep ocean crust corresponding to the intervals 
- 2.8, 3.7 - 5.9 Myr (see Table 1). They concluded that 
these ^"Fe atoms could only be produced by a SN ex- 
plosion, which they proposed took place about ~ 5 Myr 
ago at 13 ^ 30 pc, causing the excess of ^°Fe in the sec- 
ond layer. The youngest layer results were tentatively 
explained as due to a background of radioactive iron in 
the solar neighborhood. A reanalysis of these data came 
to similar conclusions!^, but attributed the presence of 
^"Fe in the younger layer to biomixing. Here we propose 
that the origin of the ^°Fe atoms are the Sco-Cen SNe. 
As we show below, both the amplitude and timing of 
their expected deposition rate are in the right range to 
explain the observed excess. 

Gas ejecta from a SN explosion cannot easily reach 
the Earth unless the pressure from the SN blast front 
is larger than the ram pressure of the solar wind at the 
Earth orbit. For an isolated SN, whose front reaches the 
Solar System driven by momentum conservation 0, this 
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TABLE L Predictions and measurements of Fe ex- 
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FIG. 1: Evolution ol the total distance between the Sun 
and Sco-Cen subgroups during the last 11 Myr. For each 
subgroup, only the epoch during which SNe were being formed 
is shown (see also Fig. 1 of Ref. |Q which shows the projection 
onto the Galactic Plane of the positions of the Earth and the 
Sco-Cen association) 



roughly defines a cut-ofF distance of 15 — 100 pc (depend- 
ing on the geometrical configuration of the SN, the Sun 
and the Earth) . This is also true in the case of multiple 
SN, though in that case the second and subsequent SNe 
explode in a very low density medium, making their ex- 
pansion follow the Sedov regime ||l[ |lC, 11 1 up to distances 
comparable with Local Bubble size ~ 200 pc. Although 
this would make it difficult for gaseous debris from all 
but the closest Sco-Cen SNe to penetrate the heliosphere 
up to the Earth's radius, observations show that most of 
the iron in the Local Bubble is condensed forming dust 
[0, [T^ . The interaction of interstellar dust with the so- 
lar wind and magnetic field is a complex problem, which 
strongly depends on the size of the dust particles [p^, 
but interstellar dust containing iron has been found at 
the Earth's orbit jl^, and it seems reasonable to assume 
that most, or at least a large fraction of the iron reaching 
the heliosphere traveling in a SN blast front would have 
reached the Earth's orbit. p^. 

Table 1 presents the typical distance and expected 
number of Sco-Cen SNe in each of the corresponding time 
intervals. The expected surface density of ^°Fe (corrected 
for in situ decay) deposited in a layer N{Al) can be es- 
timated M as: 



N{Al) ^ 4.1 X IO^A^sn/ 
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where A^sn is the number of SNe, which are assumed to 
happen at a typical distance D{Al) during the time 



interval covered by each of the sediment layers, / is the 
uptake factor that estimate as 1/100, and Meopois the 
expected ^"^Fe yield by a SN. To compare with the results 
of Ref. 0, we divide N{Al) by the At covered by each 
layer, which yields the flux 0sn (cm~^ Myr ~^)presented 
in Table 1. 

This estimation has several sources of uncertainty. The 
Fe yield can vary from 10~^ Mq to 10~® Mq depend- 
ing on the SN mass (type II SN, |l|, |l|, |l^, |l^ |§ ) . 
Also, the uptake factor /, which represents the fraction of 
^°Fe atoms present in the ocean which is deposited in the 
crust, has large and difficult to estimate uncertainties 
To compound these two error sources, the positions and 
time of the explosions can have a considerable scatter, 
and even the total number of SNe has a Poisson uncer- 
tainty of 24%. Therefore, we consider the values of (/)sn 
in Table 1 to be, at best, order-of-magnitude estimates, 
and that is reflected in the error bars assigned to our 
predictions in Fig 2. Table 1 also presents the inferred 
average fluxes (pohs measured by Ref. 0, and the ex- 
pected background level (t>h based on a measurement of 
a 13 Myr old core. Note that in Fig 2, we have added 
the ^"Fe background rate (j)h to the contribution from 
the Sco-Cen SNe 0sn , and that the error bars correspond 
only to the uncertainty in (j)sN, since the above reference 
does not provide an error estimate for the background. 

The agreement is excellent for the first and youngest 
layer, which has the highest signal-to-noise. The flux 
measured in the second layer is a factor 4 higher than 
our prediction, but there is an ample overlap between 
the error bars of both quantities, and thus the results 
can be considered consistent. Regarding the third, old- 
est layer, it is not clear whether there is any signal above 
the background, so our prediction of a very small flux 
is also compatible with the observation . It should also 
be noted that no signal of SN origin was detected in the 
same layers for another isotope ^'^Mn |Q, ||] , which is again 
consistent with our scenario since the predicted ^^Mn de- 
position rates are much lower than the background of cos- 
mogenic origin. We therefore conclude that the Sco-Cen 
SNe are enough to explain the excess of ®°Fe in the deep 
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FIG. 2: Comparison between the Fe deposition rate pre- 
dicted from the Sco-Cen SNe and the measurements of deep 
ocean crust samples. The horizontal continuous lines/empty 
squares and the shadowed boxes represent respectively the 
data and errors of Ref. , while our predictions and associ- 
ated error estimates are represented by the filled squares and 
error bars. The dashed lines correspond to the background 
estimates of Ref. M 



ocean crust. It will be very interesting to obtain crust 
data with better "temporal" resolution to identify indi- 
vidual SN explosions; including in the search other SN- 
created radioisotopes with low backgrounds and proper 
decay rates would allow to pin down the mass range of the 
progenitors and the distance at which the SNe exploded. 

Several authors |21, 2^, Q have proposed that the 
explosion of a nearby SN could have caused one or more 
of the massive extinction events in the fossil record. At 
distances larger than a few pc, the only component of the 
SN emission capable of performing serious damage to the 
biosphere is the charged cosmic ray radiation |^ (for 
instance, hard UV radiation from the SN shock break- 
out would produce flux levels above the atmosphere of 
Fuv ~ 10^(50pc/D)2 ergs cm-2 during 1 day ||^, a 
level smaller than the amount of similar radiation re- 
ceived from the Sun [^). A strong increase in the flux 
of cosmic rays reaching the upper levels of the atmosphere 
speeds up the production of NO, which catalytically de- 
stroys large amounts of ozone molecules . Assuming 
that the energy released by the SN in cosmic rays is 10^° 
ergs, the time integrated flux or fluence reaching Earth 
wiU be I2I: 



At « 2.2 X 10' 
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that the cosmic rays would travel by diffusion in a 
random Galactic magnetic field with a scale-length of 1 
pc and that, therefore, the above flux would be spread 
over a period of At w 3(D/1 pc)^ yr. However, the as- 
sumption of randomness is not valid in this case: the 
magnetic field in the outer 'shell' of the Local Bubble is 
probably coherent on large scales and very weak or inex- 
istent in its interior , being pushed out by the effects 
of the SN explosions. Therefore, the amount of cosmic 
rays reaching the Earth in the event of a nearby SN ex- 
plosion would strongly depend on the relative position 
of the SN with respect to Earth and the Local Bubble 
shell. If both the SN and the Sun are within the cavity, 
the cosmic rays will probably travel unhindered and hit 
the Earth spread over a At ~ r, where t corresponds 
to their emission period, estimated to be in the range 
T ~ 10 - 10 il, io|- If the SN were in the outside, 
the Bubble shell would serve as a "conductor" for the 
cosmic rays if the Sun were also outside the Bubble, or 
as an "insulating" shield, if the Sun were inside. 

This means that the value of the time interval At is 
highly uncertain, and the best we can do is estimate the 
maximal CR flux produced by one of the Sco-Cen SNe 
by taking the fluence above and a At similar to the lower 
limit of the typical cosmic ray acceleration time r ^ 10 
yr. Using the minimal distance for the Sco-Cen SNe D > 
40pc we flnd that (pen. < 1-4 x lO^ergs cm~^yr^ ^ , within 
a factor 2 of the flux estimated by for a 10 pc event 
using the random magnetic fleld assumption. It has been 
estimated |3^ that such an event would lead to a ozone 
depletion of 60% at high latitudes and about 20% at the 
equator. 

Therefore, the subsequent increase in the UV-B flux 
from the Sun reaching the Earth surface due to one or 
more of the Sco-Cen SNe could have caused at most a 
minor extinction |32, p3[, but would not be enough to 



provoke a major mass extinction like e.g. the Cretaceous- 
Tertiary event (see also , where the mutagenic effects 
of excess UV radiation are discussed in detail). This ex- 
tinction would particularly affect marine ecosystems, as 
flrst proposed by Ref. ; an increase on the UV-B flux 
over the ambient level can provoke a signiflcant reduction 
in phytoplankton abundance and biomass, propagating 
to at least one species of zooplankton as a secondary ef- 
fect Q. Despite the uneven distribution of the ozone 
depletion, tropical species would be more affected due 
to the higher solar angle at which they receive their UV 
dose 1 33 . Therefore, the biological signature of a SN 



It has usually been considered in the literature P4l 



explosion provoking a "clean" UV-B catastrophe would 
be a decline of ocean surface phytoplankton productiv- 
ity not associated with other causes as volcanic activity, 
climate changes or impact events. A decline in plankton 
productivity would be very difficult to detect in the fossil 
record, but it could be inferred from secondary effects, in 
particular an extinction of moUusks . 

Such an extinction affecting marine tropical, subtrop- 
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ical and temperate American bivalves characterized the 
Phocene-Pleistocene boundary, ~ 2 Myr ago |3|, |3^, |3^, 
|38t (see also Ref. Q). Two hypotheses have been ad- 
vanced to explain this phenomenon: the emergence of the 
Panama isthmus ^ and cooling due to the onset of 
Northern Hemisphere glaciations Q. However, 

this extinction episode was too rapid to be due to the 
Panama isthmus closure Q and a detailed analysis 
of the extinction patterns seems to rule out cooling as the 
]3^ , 40 . A simultaneous, although slower, episode 



cause 

of extinction affected corals |44 , which are known to be 
highly susceptible to UV-B radiation This leaves 

a SN-provoked UV-B catastrophe as a possible candi- 
date for the Pleistocene-Pliocene extinction; it should 
be noted that this epoch roughly coincides with the time 
of closest approach of LCC (see Fig. 1), during which 
the probability of nearby SN explosions would have been 
highest. In addition, Ref. Q proposed that a SN blast 
wave ionized the local interstellar medium between 2 and 
3.6 Myr ago. 

To test this hypothesis, the time and distance at which 
individual SN explosions took place should be deter- 
mined more precisely, using geological information as 
suggested above. A coincidence in time between the SN 
expected to have strongest effects on the biosphere and 
the Pleistocene-Pliocene extinction, would strongly sup- 
port the existence of a link between both events. 
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